We report a detailed DFT study of the energetic and structural properties of the spin-crossover Co(II) complex [Co(tpy)₂]²⁺ (tpy = 2,2′:6′,2′′-terpyridine) in the low-spin (LS) and the high-spin (HS) states, using several generalized gradient approximation and hybrid functionals. In either spin-state, the results obtained with the functionals are consistent with one another and in good agreement with available experimental data. Although the different functionals correctly predict the LS state as the electronic ground state of [Co(tpy)₂]²⁺, they give estimates of the HS–LS zero-point energy difference ΔE⁰_HL (tpy)  which strongly depend on the functional used. This dependency on the functional was also reported for the DFT estimates of the zero-point energy difference ΔE⁰_HL (bpy)  in the HS complex [Co(bpy)₃]²⁺ (bpy = 2,2′-bipyridine) [A. Vargas, A. Hauser and L. M. Lawson Daku, J. Chem. Theory Comput., 2009, 5, 97]. The comparison of the ΔE⁰_HL (tpy)  and ΔE⁰_HL (bpy)  estimates showed that all functionals correctly predict an increase of the zero-point energy difference upon the bpy → tpy ligand substitution, which furthermore weakly depends on the functionals, amounting to (ΔE⁰_HL)_bpy->tpy  ≈ +2670 cm^-1 . From these results and basic thermodynamic considerations, we establish that, despite their limitations, current DFT methods can be applied to the accurate determination of the spin-state energetics of complexes of a transition metal ion, or of these complexes in different environments, provided that the spin-state energetics is accurately known in one case. Thus, making use of the availability of a highly accurate ab initio estimate of the HS–LS energy difference in the complex [Co(NCH)₆]²⁺ [L. M. Lawson Daku, F. Aquilante, T. W. Robinson and A. Hauser, J. Chem. Theory Comput., 2012, 8, 4216], we obtain for [Co(tpy)₂]²⁺ and [Co(bpy)₃]²⁺best estimates of ΔE⁰_HL (bpy) ≈ -2800 cm^-1  and ΔE⁰_HL (tpy) ≈ 0 cm^-1 , in good agreement with the known magnetic behaviour of the two complexes.

Density functional theory is applied within a supramolecular approach to the study of the guest−host interactions in [Fe(bpy)₃]²⁺@Y and their influence on the structural, energetic, and ⁵⁷Fe Mössbauer spectroscopy properties of the encapsulated [Fe(bpy)₃]²⁺ complex in the low- and high-spin states. The structures of the isolated complex and the supramolecular model used for [Fe(bpy)₃]²⁺@Y were optimized in both spin-states using different generalized gradient approximation (PBE, HCTH, OLYP) and hybrid (B3LYP*, O3LYP) functionals. The results obtained are consistent with one another and show that, in either spin-state, the structure of [Fe(bpy)₃]²⁺ shrinks and distorts upon encapsulation. Still, the structural changes experienced by the complex in a given spin-state remain limited, especially in that they do not lead to a substantial variation of the ⁵⁷Fe quadrupole splitting, whose calculated values are in very good agreement with avalaible experimental data. The decomposition of the guest−host interaction energy into its electrostatic, Pauli and orbital contributions shows that the bonding between the complex and the supercage is more electrostatic than covalent. The ability of modern functionals to accurately describe the interactions explains the remarkable consistency of the results obtained with the various functionals. In particular, although the functionals perform very differently for the determination of the high-spin/low-spin energy difference ΔE_HL^el in [Fe(bpy)₃]²⁺ and [Fe(bpy)₃]²⁺@Y, they consistently predict that the encapsulation entails a destabilization of the high-spin state with regard to the low-spin state of Δ(ΔE_HL^el) = 2500 cm^−1. Using for [Fe(bpy)₃]²⁺ the CASPT2 value of ΔE_HL^el = 3700 cm^−1 [Pierloot, K.; Vancoillie, S. J. Chem. Phys.2006, 125, 124303; Pierloot, K.; Vancoillie, S. J. Chem. Phys.2008, 128, 034104], we obtain for the high-spin/low-spin energy difference in [Fe(bpy)₃]²⁺@Y, a best ab initio estimate of ΔE_HL^el = 6200 cm^−1.

The trigonal planar geometry of the nitrogen atom in commonly used phosphoramidite ligands is not in line with the traditional valence shell electron pair repulsion (VSEPR) model. In this work, the effects governing nitrogen configuration in several substituted aminophosphines, A₂PNB₂ (A or B = H, F, Cl, Br, Me, OMe, BINOP), are examined using modern computational analytic tools. The electron delocalization descriptions provided by both electron localization function (ELF) and block localized wavefunction analysis support the proposed relationships between conformation and negative hyperconjugative interactions. In the parent H₂PNH₂, the pyramidal nitrogen configuration results from nitrogen lone pair electron donation into the σ* P — H orbital. While enhanced effects are seen for F₂PNMe₂, placing highly electronegative fluorine substituents on nitrogen (i.e., Me₂PNF₂) eliminates delocalization of the nitrogen lone pair. Understanding and quantifying these effects can lead to greater flexibility in designing new catalysts.

Whereas there are hundreds of known iron(II) spin-crossover compounds, only a handful of cobalt(II) spin-crossover compounds have been discovered to date, and hardly an in depth study on any of them exists. This review begins with an introduction into the theoretical aspects to be considered when discussing spin-crossover compounds in general and cobalt(II) systems in particular. It is followed by case studies on [Co(bpy)₃]²⁺ and [Co(terpy)₂]²⁺ (bpy = 2,2′-bipyridine, terpy = 2,2′:6′,2″-terpyridine) presenting and discussing results from magnetic susceptibility measurements, X-ray crystallography, optical spectroscopy, and EPR spectroscopy.

State-of-the-art generalized gradient approximation (GGA) (PBE, OPBE, RPBE, OLYP, and HCTH), meta-GGA (VSXC and TPSS), and hybrid (B3LYP, B3LYP*, O3LYP, and PBE0) functionals are compared for the determination of the structure and the energetics of the D₃ [Co(bpy)₃]²⁺ complex in the ⁴A₂ and ⁴E trigonal components of the high-spin ⁴T_1g( t⁵_2g e²_g ) state and in the low-spin ²E state of octahedral ²E_g( t⁶_2g e¹_g) parentage. Their comparison extends also to the investigation of the Jahn−Teller instability of the ²E state through the characterization of the extrema of C₂ symmetry of this spin state's potential energy surface. The results obtained for [Co(bpy)₃]²⁺ in either spin manifold are very consistent among the functionals used and are in good agreement with available experimental data. The functionals, however, perform very differently with respect to the spin-state energetics because the calculated values of the high-spin/low-spin energy difference ΔE^el_HL vary between −3212 and 3919 cm^-1. Semilocal functionals tend to give too large ΔE^el_HL values and thus fail to correctly predict the high-spin state as the ground state of the isolated complex, while hybrid functionals tend to overestimate the stability of the high-spin state with respect to the low-spin state. Reliable results are, however, obtained with the OLYP, HCTH, B3LYP*, and O3LYP functionals which perform best for the description of the isolated complex. The optical properties of [Co(bpy)₃]²⁺ in the two spin states are also analyzed on the basis of electronic excitation calculations performed within time-dependent density functional response theory. The calculated absorption and circular dichroism spectra agree well with experimental results.

The high-spin → low-spin relaxation in spin-crossover compounds can be described as non-adiabatic multi-phonon process in the strong coupling limit, in which the low-temperature tunnelling rate increases exponentially with the zero-point energy difference between the two states. Based on the hypothesis that the experimental bond length difference between the high-spin and the low-spin state of ~0.2 Å is also valid for low-spin iron(II) complexes, extrapolation of the single configurational coordinate model allows an estimate of the zero-point energy difference for low-spin complexes from kinetic data. DFT calculations on low-spin [Fe(bpy)₃]²⁺ support the structural assumption. However, for low-spin [Fe(terpy)₂]²⁺ the relaxation rate constant shows an anomalous behaviour in so far as it is more in line with spin-crossover systems. This is attributed to very anisotropic bond length changes associated with the spin state change, and the subsequent breakdown of the single mode model.

In the iron(II) low-spin complex [Fe(bpy)₃]²⁺, the zero-point energy difference between the ⁵T_2g(t⁴_2ge²_g) high-spin and the ¹A_1g(t⁶_2g) low-spin states, ΔE⁰_HL, is estimated to lie in the range of 2500-5000 cm^-1. This estimate is based on the low-temperature dynamics of the high-spin→low-spin relaxation following the light-induced population of the high-spin state and on the assumption that the bond-length difference between the two states Δr_HL is equal to the average value of ≈0.2 Ã…, as found experimentally for the spin-crossover system. Calculations based on density functional theory (DFT) validate the structural assumption insofar as the low-spin-state optimised geometries are found to be in very good agreement with the experimental X-ray structure of the complex and the predicted high-spin geometries are all very close to one another for a whole series of common GGA (PB86, PW91, PBE, RPBE) and hybrid (B3LYP, B3LYP*, PBE1PBE) functionals. This confirmation of the structural assumption underlying the estimation of ΔE⁰_HL from experimental relaxation rate constants permits us to use this value to assess the ability of the density functionals for the calculation of the energy difference between the HS and LS states. Since the different functionals give values from -1000 to 12000 cm^-1, the comparison of the calculated values with the experimental estimate thus provides a stringent criterion for the performance of a given functional. Based on this comparison the RPBE and B3LYP* functionals give the best agreement with experiment.

Molecular modelling and site-directed mutagenesis were used to identify eleven amino acid residues which may be involved in antagonist binding of the human tachykinin NK1 receptor. Recombinant receptors were expressed in mammalian cells using the Semliki Forest virus system. Wild type and mutant receptors showed similar expression levels in BHK and CHO cells, verified by metabolic labelling. Binding affinities were determined for a variety of tachykinin NK1 receptor antagonists in SFV-infected CHO cells. The binding affinity for GR203040, CP 99,994 and CP 96,345 was significantly reduced by mutant Q165A. The mutant F268A significantly reduced the affinity for GR203040 and CP 99,994 and the mutant H197A had reduced affinity for CP 96,345. All antagonists seemed to bind in a similar region of the receptor, but do not all rely on the same binding site interactions. Functional coupling to G-proteins was assayed by intracellular Ca2+ release in SFV-infected CHO cells. The wild type receptor and all mutants except A162L and F268A responded to substance P stimulation.

The 1.3 S biotinylatable subunit of Proprionibacterium shermanii transcarboxylase complex was fused to the C-terminus of the human neurokinin 1 receptor gene and introduced into the Semliki Forest virus expression vector pSFV1. RNA transcribed from pSFV1-NK1-biot and pSFV-Helper2 was coelectroporated into BHK cells permitting in vivo packaging of recombinant virus. Infection of BHK and CHO cells with SFV-NK1-biot virus yielded high level of the fusion receptor as detected by metabolic labeling, immunoblotting with streptavidin alkaline phosphatase and binding to substance P. Like native receptor, the biotinylated receptor fusion was able to stimulate Ca2+ mobilization in infected CHO cells,indicating functional coupling to guanine-nucleotide-binding proteins.